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We examine finite-temperature phase transitions in the two-orbital Hubbard model with different 
bandwidths by means of the dynamical mean-field theory combined with the continuous-time quan¬ 
tum Monte Carlo method. It is found that there emerges a peculiar slope-reversed first-order Mott 
transition between the orbital-selective Mott phase and the Mott insulator phase in the presence 
of Ising-type Hund’s coupling. The origin of the slope-reversed phase transition is clarified by the 
analysis of the temperature dependence of the energy density. It turns out that the increase of 
Hund’s coupling lowers the critical temperature of the slope-reversed Mott transition. Beyond a 
certain critical value of Hund’s coupling the first-order transition turns into a finite-temperature 
crossover. We also reveal that the orbital-selective Mott phase exhibits frozen local moments in the 
wide orbital, which is demonstrated by the spin-spin correlation functions. 
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I. INTRODUCTION 

Coexistence of strongly and weakly correlated elec¬ 
trons has been one of the intriguing subjects in condensed 
matter physics. Materials in which more than one orbital 
is active near the Fermi level have exhibited interesting 
properties and their main origin is believed to be the 
coexistence of electrons with different degrees of correla- 
tionsi. In the multiorbital system, correlations between 
electrons and Hund’s coupling have been known to show 
rich phenomena in the presence of the orbital degree of 
freedom. In case that the degeneracy between active or¬ 
bitals is lifted by the difference of their bandwidths^rJ^ 
or crystal-field splittingi^^— ^ the degree of effective cor¬ 
relations in each orbital becomes different. One promi¬ 
nent consequence of different degrees of correlations is the 
orbital-selective Mott phase (OSMP), where electrons in 
some orbitals are totally localized due to the Mott physics 
while other orbitals are still occupied by itinerant elec- 
tron a^^d^ . Here Hund’s coupling tends to intensify the 
difference between orbital o^d° . 

The coexistence of strongly and weakly correlated elec¬ 
trons is also believed to play an important role in two- 
dimensional materials including strong spatial fluctua- 
tions^iiS^. In such a system spatial correlations and cor¬ 
responding momentum-space anisotropy of correlations 
are the key elements to host the coexistence^. Thanks to 
the recent numerical developments in the cluster dynami¬ 
cal mean-field theory (DMFT)2i“— , it is known that spa¬ 
tial fluctuations modify qualitatively finite-temperature 
behaviors of the correlation-driven metal-insulator tran¬ 
sitions; this has been revealed by the comparison with the 
single-site DMFT— neglecting spatial fluctuations. Spa¬ 
tial correlations turn out to reduce greatly the ground- 
state entropy of the paramagnetic Mott insulator (MI) 
at low temperatures and accordingly, the itinerant bad 
metallic phase dominates in the region of relatively high 


temperatures near the transitioii^^. 

It is natural to anticipate such prominent changes in 
finite-temperature transitions for multiorbital systems. 
In spite of extensive studiesi2rd^i^“i2^, the temperature 
dependence of the transitions in the two-orbital Hub¬ 
bard model still lacks a thorough understanding. The 
principal purpose of this work is to investigate the finite- 
temperature nature of the transitions in two-orbital sys¬ 
tems with emphasis on the effects of Hund’s coupling. 

In this paper we investigate the two-orbital Hubbard 
model by the DMFT combined with the continuous-time 
quantum Monte Carlo (CTQMC) method. In the model 
we find the slope-reversed Mott transition in the pres¬ 
ence of Ising-type Hund’s coupling for two orbitals of 
different bandwidths. We also observe that the drastic 
changes in the phase transition between the OSMP and 
the MI phase are induced by the variation of the Hund’s 
coupling strength. The analysis of the hysteresis behav¬ 
ior of local magnetic moments determines the location 
of the critical end points, which reveals that the critical 
temperature tends to reduce as the Hund’s coupling is 
increased. Eventually the hysteretic behavior disappears 
at a certain value of Hund’s coupling and the system ex¬ 
hibits only a crossover between the OSMP and the MI 
phase. We also compute the spin-spin correlation func¬ 
tion for both orbitals and find the formation of the local 
frozen moments for itinerant wide-orbital electrons in the 
OSMP. 

This paper is organized as follows: In Sec. H we give a 
brief description of the two-orbital Hubbard model and 
the numerical method. Section HI is devoted to the pre¬ 
sentation of the numerical results, which include finite- 
temperature phase diagrams, spectral functions, hystere¬ 
sis of local magnetic moments, energy densities, effects 
of Hund’s coupling, and spin-spin correlation functions. 
The results are summarized in Sec. IV. 
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II. MODEL AND METHODS 

We consider the Hamiltonian 
'^ = - + h.C.) - 

{ij)a.G iota 

+ - J5acr-)nii^na„', ( 1 ) 

ioi i<7<7' 

for two orbitals a = 1 and 2. Here, c^^ai^laa) i® 
nihilation(creation) operator of an electron with spin a 
at site i and orbital a. In each orbital electrons move on 
the infinite-dimensional Bethe lattice corresponding to 
a noninteracting semi-circular density of states (DOS), 
p°(w) = {2/TTDa)y/l — {ijj/DaY, with the half band¬ 
width Da = 2ta and interact with each other via the 
intra- and inter-orbital Coulomb interactions U and U' 
and Hund’s coupling J. We investigate the half-filled 
system with chemical potential fx = 317/2 — 5J/2, and 
also choose D 2 = 2Di and U' = U — 2J. The half band¬ 
width Di of the narrow orbital is taken as an energy unit 
throughout this paper. Here we disregard spin-flip and 
pair-hopping terms, which is appropriate for the study of 
the anisotropic Hund’s coupling model. Our model serves 
as a natural generalization of the Ising-spin Kondo lat¬ 
tice model which is useful in interpreting experimental 
results for various materials, e.g., pyrochlore oxides^i^ 
and URu 2 Si 2 ^. Such physics is understood mainly in 
terms of the large-anisotropy effects on the localized mo¬ 
ments. 

We employ the DMFT combined with the CTQMC 
method through the hybridization expansion algo- 
rith m^"^d°d^ . Typically, the statistical sampling of 10® 
Monte Carlo steps is performed, which turns out to be 
sufficient for statistically reliable numerical results. 

III. RESULTS 

A. Finite-Temperature Phase Diagram 

The main result of this work is the emergence of a 
slope-reversed Mott transition accompanied by the dras¬ 
tic change in the behavior of finite-temperature phase 
transitions, which is driven by the variation in Hund’s 
coupling. Figures [I] shows phase diagrams on the temper¬ 
ature versus interaction strength plane for various Hund’s 
coupling strengths. In the presence of orbital degrees of 
freedom, generally, we have two successive phase transi¬ 
tions, one from the Fermi-liquid (FL) phase to the OSMP 
and the other from the OSMP to the MI phase. The 
transition between the FL phase and the OSMP inher¬ 
its the shape and energy scale of the coexistence region 
in the single-orbital model. In Fig. [T](a) and (b), on the 
other hand, the coexistence region of the OSMP-to-MI 
phase transition is quite interesting. First of all, the 
slope of the phase-transition line is opposite to that in 
the single-orbital case shown in the inset of Fig. [IJa). 


The slope-reversed Mott transition was reported in the 
two-dimensional systems and its origin was attributed to 
spatial modulations^^. Here it is noted that our system 
is an infinite-dimensional one without any spatial fluc¬ 
tuations. We can also find that the critical temperature 
associated with the slope-reversed transition is consider¬ 
ably enhanced. 

The effects of Hund’s coupling are rather drastic on 
the slope-reversed Mott transition. When we increase the 
Hund’s coupling strength, the slope-reversed Mott tran¬ 
sition becomes a finite-temperature crossover, implying 
a continuous transition at zero temperature. A similar 
change in the zero-temperature transition was reported 
in an effective low-energy model^; our result reveals that 
the zero-temperature result reflects the change from the 
slope-reversed transition to a crossover at finite temper¬ 
atures. In addition, the region of the OSMP, which is 
present between the two transitions, becomes wider for 
larger Hund’s coupling strength, from which we can infer 
that Hund’s coupling plays the role of a ‘band decou¬ 
pler’—. It is also found that for very small Hund’s cou¬ 
pling strength, J = C//32, the coexistence regions of the 
two transitions overlap significantly with each other. In 
Fig. [T](d) we also plot the existing results obtained from 
exact diagonalization (ED}^d^ and Hirsch-Fye quantum 
Monte Carlo (HF-QMC) ^°’®® , which are reasonably con¬ 
sistent with our numerical results. 

The reversed slope of the phase-transition line is a dis¬ 
tinctive feature. In a conventional Mott transition the 
localized MI phase dominates the itinerant FL phase in 
the region of high temperatures near the phase transi¬ 
tion; this is mainly due to the extensive entropic con¬ 
tribution of the MI phase compared with the very small 
ground-state entropy in the FL phase. Similarly to the 
slope-reversed transition in two-dimensional systems, the 
origin of which is the significant entropy reduction of the 
MI phase by the short-range correlations^, the slope- 
reversed Mott transition in the two-orbital system can 
be understood in terms of the entropy of the MI phase: 
It is expected to reduce considerably through ferromag¬ 
netic correlations between electrons in different orbitals 
by Hund’s coupling. Another important aspect in the 
two-orbital system is that instead of the FL phase, the 
OSMP competes with the MI phase near the transition. 
The OSMP, in which electrons are partly localized, has 
higher entropy than the FL, and accordingly it is more 
likely to dominate the MI phase at high temperatures 
to yield the reversed slope of the transition line. We 
will give a detailed analysis in a later subsection, where 
the temperature-dependence of the energy density is dis¬ 
cussed. 


B. Spectral Function and Self-Energy 

The local spectral function of each orbital, which 
can be evaluated via an analytic continuation to the 
real-frequency domain by the maximum entropy method 
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FIG. 1. (Color online) Phase diagrams on the plane of temperature T and interaction strength U for (a) J = !7/32, (b) 
J = 17/16, (c) J = U/8, and (d) J = t//4. In (a) and (b), open and filled symbols correspond to lower and upper transition 
interaction strengths, respectively. The transitions for narrow and wide orbitals are denoted by (blue) squares and (red) circles, 
respectively. The inset in (a) exhibits transition interaction strengths of the single-orbital Hubbard model. In (c) and (d) 
(blue) open squares and (blue) filled squares represent lower and upper transition interaction strengths, respectively, of the 
narrow orbital transition and (red) filled circles the phase boundary within the crossover region between the OSMP and the 
MI phase. The (black) pluses J+] as well as the dashed line represent the result of finite-temperature ED (Ref. [l^ while that 
of zero-temperature ED (Ref. |y) is represented by (black) crosses [x]. The results from the HE-QMC method in Refs. [1(1 and 
1^ are marked by a (black) open triangle and asterisks [*], respectively. The diamonds represent the critical end points of the 
corresponding first-order transitions. 


(MEM), characterizes conveniently the feature of each 
phase in the phase diagram. In the left panel of Fig. [21 
the spectral functions of three different phases are shown 
for J = 17/16. In the FL phase with U = 3.0, the spectral 
function exhibits clearly a coherent peak, which satisfies 
the Luttinger theorem. On the other hand, the coher¬ 
ent peak disappears and Mott gaps develop for both or¬ 
bitals in the MI phase. For the intermediate interaction 
strength corresponding to OSMP, the narrow orbital is 
gapped while the wide one still remains itinerant. It is 
remarkable that the spectral function of the wide orbital 
deviates substantially from the noninteracting DOS at 
the Fermi level. The violation of the Luttinger theo¬ 
rem implies the finite lifetime of wide-orbital electrons 
at the Fermi level. The finite-scattering amplitude of 
the wide-orbital electron at the Fermi level can be veri¬ 
fied by the finite offset in the imaginary-part of the self¬ 


energy, as shown in Fig.|2{e). Similar evidences were also 
reported for the non-Fermi-liquid nature of the OSMP 
which crosses over to the MI phas o^^i^^ . 

C. Local Magnetic Moments 

The first-order transition between the OSMP and the 
MI phase is demonstrated by the hysteresis behavior of 
physical quantities such as the local magnetic moment. 
In Fig. Ha) we plot the local magnetic moment of the 
wide orbital as a function of U for different tempera¬ 
tures. As the interaction strength is increased, electrons 
become more localized and the average local moment in¬ 
creases monotonically. Over a finite region of the in¬ 
teraction strength we can observe the hysteresis of the 
local spin magnetic moment, which implies the coexis- 
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FIG. 2. (Color online) Spectral functions calculated via the 
maximum entropy method [(a) to (c)] and imaginary-part self¬ 
energies of wide-orbital electrons [(d) to (f)] at T = 1/200 for 
J/U = 1/16. From top to bottom, the interaction strength 
corresponds to U — 3.00 [(a) and (d)], 3.60 [(b) and (e)[, and 
4.00 [(c) and (f)[. In (a) to (c), the (blue) solid and the (red) 
dashed lines represent the spectral functions of narrow and 
wide orbitals, respectively. For comparison, the noninteract¬ 
ing density of states, marked with the (black) dot-dashed line, 
is also shown in (a) both orbitals and (b) wide orbital. 

fence of the two phases. As shown in Fig. [31 we can es¬ 
timate two transition interaction strengths and U ^2 
from the minimum and the maximum values of U, respec¬ 
tively, showing the coexistence. The coexistence region 
shifts to the stronger interaction region with the increase 
of the temperature, resulting in the reversed slope of the 
phase-transition line. 

Using the above hysteresis, we can also estimate the 
position of the critical end point of the slope-reversed 
Mott transition. From the numerical data, we obtain the 
maximum difference of the local moments for the two 
solutions (MI phase and OSMP) in the coexistence region 

HSz2) = Max [(S'^2 )mi - (*S'^2)osmp] (2) 

at each temperature. In the plot of as a func¬ 

tion of T, the T-axis cut gives the critical temperature 
T>, as shown in Fig. |31[b). The hysteresis data provide 
the interaction strength Umax, where S{S‘^ 2 ) reaches the 
maximum, and the extrapolated value of Umax to T = 
gives the interaction strength of the critical end point. 

[See Fig. |31[c).] We have thus determined the location 
of the critical end points for both first-order transitions, 
which are plotted in Fig. |TJ 


FIG. 3. (Golor online) (a) Squared magnetic moments 
of wide-orbital electron as a function of the interaction 
strength at various temperatures for J/U = 1/16. From 
top to bottom, corresponding temperatures are T = 
1/33,1/40,1/50, 1/67, 1/100, and 1/200. For better compari¬ 
son, the data for T — 1/40,1/50, 1/67, 1/100, and 1/200 are 
shifted downward by 0.02,0.04,0.06,0.08, and 0.10, respec¬ 
tively. Lines are guides to the eye. (b) Maximum differences 
of the squared magnetic moments between two solutions in 
the coexistence region at given temperatures. The solid line 
corresponds to the least-square fit of the data. The critical 
temperature T/^ is estimated by the T-axis cut of the extrap¬ 
olated line, which is denoted by the (red) dashed arrow, (c) 
Interaction strength Umax at which the difference reaches the 
maximum. The critical interaction strength is estimated 
by the the extrapolation of the least-square fit [(blue) solid 
line[ to the critical temperature [(red) vertical dashed arrow). 


D. Origin of Slope-Reversed Mott Transitions 

The investigation of the temperature-dependence of 
the total energy density sheds light on the origin of the 
slope-reversed transitions between the OSMP and MI 
phase. Figure 0] represents the total energy density e 
as a function of temperature for the three phases, where 
e is defined to be 



with N being the number of lattice sites. 

In the localized MI phase the total energy density is 
nearly constant, which reflects the fact that the entropy 
is insensitive to the temperature at low temperatures. On 
the other hand, the itinerant FL phase gives a monotonic 
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FIG. 4. (Color online) Total energy density as a function of 
temperature with J/U = 1/16 for three different interaction 
strengths U = (a) 3.0, (b) 3.6, and (c) 5.0, which correspond 
to the FL phase, the OSMP, and the MI phase, respectively. 
The inset in (a) exhibits behavior of the Fermi-liquid phase 
at low temperatures. 
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increase in the total energy density with the increase of 
the temperature. As shown in the inset, the increas¬ 
ing behavior is consistent with behavior at low tem¬ 
peratures. Interestingly, in the OSMP the total energy 
density also increases as the temperature is increased 
as in the FL phase. Such an increase in eoSMP makes 
the OSMP more favorable compared with the MI phase 
through the additional contribution to the entropy at fi¬ 
nite temperatures. Here it is noted that the temperature 
dependence in the OSMP shows the superlinear behavior, 
eoSMp('F) « AT'^ with 7 > 1. 

We can also see that the residual entropy of the ground 
state is In 2 per site in the OSMP and the MI phase. 
(Note that the Boltzmann constant has been absorbed 
in the temperature T.) In the OSMP, only the electrons 
in narrow orbitals are localized and the degree of free¬ 
dom for their spins gives the residual entropy In 2. In 
the MI phase, on the other hand, electrons in both nar¬ 
row and wide orbitals are localized. Nevertheless, the 
Ising-type Hund’s coupling makes the ground state of 
the local Hamiltonian be still two-fold degenerate, com¬ 
posed of I t; t) and I J,), where \a] (3) describes the state 
with a spin-a electron in the narrow orbital and a spin-,5 
electron in the wide orbital. 

Based on these results, we can construct a generic 
phase boundary between the OSMP and the MI phase. 
Suppose a zero-temperature quantum phase transition 



FIG. 5. (Golor online) (a) Schematic representation of the 
ground-state energy densities of the OSMP and the MI phase, 
which are represented by the (blue) solid and the (red) dashed 
lines, respectively. denotes the transition interaction 

strength at zero temperature, (b) Variations of the free- 
energy densities of the OSMP [(blue) solid line] and the MI 
phase [(red) dashed line] for U > as the temperature in¬ 
creases. Phase transition from the MI phase to the OSMP oc¬ 
curs at Tp only for U > U^. The inset displays the resulting 
transition temperature as a function of U, which reproduces 
well the slope-reversed Mott phase transition. 


between the OSMP and the MI phase takes place at 
U = U^Q, where the ground state energies of two phases 
cross, as illustrated in Fig. [51[a). We use the relation 


s{T) 


s{T = 0) -b 


r'^ dT' de{T’) 

Jo dT' 


(4) 


to estimate the entropy density s{T) at low tempera¬ 
tures. The resulting temperature dependencies of the 
free-energy density for the two phases are given in the 
form 

/oSMP = Eg.OSMP — Tln2-, (5) 

7-1 

/mi = £g,Mi — Tin2 . (6) 

Both phases have the same residual entropies while the 
OSMP has additional free-energy gain, shown in the third 
term of Eq. ©• This contribution originates from the su¬ 
perlinear temperature dependence of the energy density 
and the corresponding entropy gain in the OSMP at fi¬ 
nite temperatures. For U < C//^, /osmp is always lower 
than /mi since £g,oSMP < £g,Mi- For U > on the 
other hand, £g,oSMP > £g,Mi and there occurs a phase 
transition at the temperature given by 


T> 

p 


7-1 

A7 


(£g,OSMP 



(7) 


Below Tp the MI phase has lower free energy while the 
increase of temperature above Tp induces a transition to 
the OSMP phase. Near the zero-temperature transition 
interaction strength U^, the ground-state energy differ¬ 
ence £g,Mi — £g,oSMP is expected to be linearly propor¬ 
tional to ?7 — U^, resulting in the following dependence 
of the transition temperature Tp ^ {U — The 

inset of Fig. [5l[b) represents a generic phase transition 
line between the OSMP and the MI phase, which turns 
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FIG. 6. (Color online) Phase diagram on the plane of the 
interaction strength U and the Hund’s coupling strength J/U 
at the temperature T — 1/200. Filled and open (blue) squares 
indicate upper (17/5) and lower {U^) transition interaction 
strengths, respectively, of the narrow-orbital transition. For 
J/U < 0.1, the (red) Filed and open circles represent upper 
(17/5) and lower (!7/^) transition interaction strengths of the 
wide-orbital transition. For J/U > 0.1, the crossover points 
between the OSMP and the MI phase are marked by the (red) 
filled circles. Lines are guides to the eye. The transition 
interaction strengths for the narrow and the wide orbitals, 
which are estimated in the limit J ^ ta, are represented 
by (green) dotted and (blue) dash-dotted lines, respectively. 
The transition interaction strength for J = 0 is also marked 
by a (black) asterisk. Inset: Critical temperature of the wide- 
orbital transition versus J/U. 

out to be slope-reversed. The resulting phase transition 
line also reproduces well the sublinear dependence of Tp 
on U — U^, which is observed in Figs.dja) and (b). 

E. Effects of Hund’s Coupling 

In Fig. [6] we summarize the effects of Hund’s coupling 
on the transitions by plotting various transition interac¬ 
tion strengths versus J/U at T = 1/200, which is the 
lowest temperature considered. For J/U = 1/64, the 
system appears to undergo a single transition without 
the OSMP. For larger values of J/U, we can observe two 
separate transitions, and the region of OSMP expands 
gradually with the increase of J/U. It is also notable that 
the critical interaction strengths associated with both or¬ 
bitals tend to decrease as Hund’s coupling grows. 

Following the Hubbard criterion for the Mott transi¬ 
tion, which is extended for the multi-orbital modelaii^, 
we can simply estimate critical interaction strength. In 
the extremely localized atomic limit {ta = 0), the charge 


excitation gap is given by 

Aatom = [Eg{N + 1) - Ep{N)] - [Ep{N) - Eg{N - 1)] 
= (1 + J/U)U , (8) 

where Eg(n) is the ground-state energy with n electrons. 
The gap is reduced by the introduction of the kinetic en¬ 
ergy, and at the critical interaction strength the reduced 
gap vanishes: 

0 = Aatom - , (9) 

where W is the estimate of the average kinetic energy. 
For J ta the charge excitations in both orbitals are 
hybridized with each other. Accordingly, the single Mott 
transition arises in this limit. The charge excitations in 
both orbitals make contribution to the kinetic energy, 
yielding the estimate 2y/Dl -b for the average kinetic 
energy; this results in the enhanced critical interaction 
strength. In the opposite limit J ^ ta, in contrast, 
orbital fluctuations are strongly suppressed and charge 
excitations in the two orbitals are not hybridized with 
each other. The average kinetic energy of the orbital a 
reduces to the bare bandwidth 2D a, leading to the two 
transition interaction strengths 


u<= 

“ 1 + J/U ' 

(10) 

^ 2 H 2 

(11) 

{7> = 

" l + J/U' 


which generally decrease with J/U. 

The above estimates of the transition interaction 
strengths are qualitatively consistent with our numeri¬ 
cal data. The Hund’s coupling decouples the excitations 
in two different orbitals, and the transition interaction 
strengths of the two orbitals begin to be separated as the 
Hund’s coupling strength is raised. The corresponding 
OSMP region becomes enlarged in the phase diagram. 
Hund’s coupling thus plays the role of ‘band decoupler’. 
The interpolation between J = 0 and the limit ot J ^ ta 
clearly shows that the transition interaction strength is 
a decreasing function of J/U. This is a characteristic of 
the half-hlled system and different behaviors in general 
fillings were reported in several work3^ i^°d^~ — . 

The inset of Fig.IHlshows that the critical temperature 
T/J of the wide-orbital first-order transition reduces as 
J/U is increased. Above a certain value of J/U, which 
turns out to be between 1/16 and 1/8, we cannot find 
the transition down to T = 1/200, the lowest temper¬ 
ature considered, only to observe crossover phenomena. 
We presume that the critical temperature of the Mott 
transition continues to diminish as J/U increases and 
eventually becomes zero between J/U = 1/16 and 1/8; 
this explains the drastic change in transition nature from 
the first-order to crossover. 
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FIG. 7. (Color online) Imaginary-time dependence of the 
spin-spin correlation function at T = 1/200 with J = [7/16 
for (a) the narrow and (b) the wide orbitals. From bottom 
to top, the interaction strength is given by [/ = 3.1, 3.3, 3.6, 
3.8, and 4.0 . 


F. Spin-Spin Correlation Function 

We next investigate the spin-spin correlation function 

4“)(r)^(5:(r)5:(0)) (12) 

for orbital a = 1,2. The spin-spin correlation function 
can give a signal for the formation of the frozen local 
magnetic moment. The long-term memory in the corre¬ 
lation function is proportional to the magnitude of frozen 
moments. Figure [7] represents the spin-spin correlation 
function of the narrow and the wide orbitals at T = 1/200 
for J = [7/16 and various interaction strengths. 

In the FL phase, (r) for both orbitals shows 1 /r^ 
scaling for imaginary time r sufficiently far from both 
0 and /3. In the OSMP, however, we can find the for¬ 
mation of the frozen local moment in the itinerant wide 
orbital (a = 2), which exhibits the long-term memory in 
C'gg(T). (See the data for U = 3.6.) In comparison with 
the moment of the narrow orbital, that of the wide or¬ 
bital is not fully developed in magnitude. Via the second 
transition, the frozen moment of wide orbital is fully de¬ 
veloped as well, and the system enters the MI phase. In 
the OSMP, we presume that not only the local moment 
of the narrow orbital but also the frozen moment of the 
wide orbital can enhance the scattering amplitude of itin¬ 
erant electrons in the wide orbital, which is observed in 
Fig. He). This itinerant phase in the wide orbital is a 
simple example of ‘frozen-moment’ metal at half filling. 


FIG. 8. (Golor online) Local spin susceptibility as a function 
of interaction strength at T = 1/200 with J = 17/16. The 
data for the narrow and the wide orbitals are marked by (blue) 
squares and (red) circles, respectively. 

Similar phases were observed at other fillingsi^i^. 

The local spin susceptibility, defined to be 

xL1= rrfr(55(r)5^(0)), (13) 

Jo 

is also shown in Fig. H Two successive first-order transi¬ 
tions are clearly observed. The intermediate OSMP has 
moderate values of X\oc i which provides another signature 
of frozen local moment. Such two-stage saturation of the 
local susceptibility was reported earlier and Hund’s cou¬ 
pling was also emphasized as an origin of the formation 
of the local moments in itinerant components^. 

IV. CONCLUSION 

We have found the slope-reversed Mott transition in 
the two-orbital Hubbard model with Ising-type Hund’s 
coupling, in which two orbitals have different band- 
widths. The reversed slope of the phase-transition line 
between the OSMP and the MI phase can be understood 
in terms of entropy contributions which are closely re¬ 
lated to the anisotropy in the Hund’s coupling. The 
analysis of the temperature dependence of the energy 
densities together with the residual entropy has given a 
successful explanation of a generic slope-reversed transi¬ 
tion between the OSMP and the MI phase. We have 
also observed drastic changes in transition nature be¬ 
tween the OSMP and the MI phase as the Hund’s cou¬ 
pling strength is varied. As the Hund’s coupling strength 
increases, the first-order transition turns into a hnite- 
temperature crossover, implying a quantum phase tran¬ 
sition at zero temperature. Such a drastic change in the 
transition nature is apparently induced by the diminish¬ 
ing critical temperature of the first-order transition be¬ 
tween the OSMP and the MI phase. Finally, the frozen 








local moments have been observed for the wide orbital in 
the OSMP. 
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